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Abstract 

We report the first measurement of the top quark mass using the decay length technique in pp 
collisions at a center-of-mass energy of 1.96 TeV. This technique uses the measured flight distance 
of the b hadron to infer the mass of the top quark in lepton plus jets events with missing transverse 
energy. It relies solely on tracking and avoids the jet energy scale uncertainty that is common to all 
other methods used so far. We apply our novel method to a 695 pb _1 data sample recorded by the 
CDF II detector at Fermilab and extract a measurement of mt = 180.7 ^13 4 (stat.) ± 8.6 (syst.) 
GeV/c 2 . While the uncertainty of this result is larger than that of other measurements, the domi- 
nant uncertainties in the decay length technique are uncorrelated with those in other methods. This 
result can help reduce the overall uncertainty when combined with other existing measurements of 
the top quark mass. 

PACS numbers: 14.65.Ha,12.15.Ff 
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I. INTRODUCTION 



A precise determination of the top quark mass (m t ) is an important goal of high-energy 
physics. The uncertainty on m t is a dominant uncertainty in global standard model (SM) 
fits for the mass of the unobserved Higgs boson. A precision measurement of m t constrains 
the allowed Higgs mass values within the SM. It will tell us where to look for the Higgs 
and test whether it is the predicted SM particle or not after a signal has been established. 
Recently, significant progress has been made in reducing the uncertainty in measurements of 
™>t [ll- Unfortunately, the most precise of the currently employed techniques are all limited 
by the same systematic uncertainty, the calorimeter jet energy scale. 

Some of the authors of this paper have developed a novel method to measure m t using 
the transverse decay length of b hadrons from top decays [2j. The method exploits the fact 
that top auarks produced in pp co.lisions at = 1.96 TeV are produced uear.y at rest aud 
decay almost instantaneously [3|] to a relatively light bottom quark and a much heavier W 
boson. In the rest frame of the top quark, the relativistic boost given to the bottom quark 
as a consequence of the top quark decay can be written simply as follows: 

= mt + m 2 b -m 2 w ^ Q ^ 
2m t m b m b 

where % = [1 — (f^/c 2 )] -1 / 2 and the approximation makes use of the fact that m t 3> mj. 
The mass of the top quark, therefore, is strongly correlated with the velocity imparted to 
the b quark and the subsequent b hadron after fragmentation. Thus, the average momenta 
of the b hadrons from top decays can be used to infer the mass of the top quark. In this anal- 
ysis, rather than measuring the average momenta, we simply measure the highly-correlated 
average transverse decay length of the b hadrons, which we denote (L 2 d)- Furthermore, we 
do not analytically solve for m t from (L 2 d), but as detailed in Section fVIII we establish the 
functional dependence of m t on (L 2 d) using Monte Carlo (MC) simulations of signal and 
background events. 

This technique relies on track reconstruction to determine precisely the decay length. The 
calorimeter information is used only for the selection of event candidates. Consequently, the 
uncertainty on the measurement due to the jet energy scale is negligible. In this paper, we 
present the first measurement of m t using the decay length technique. We apply this new 
method to pp collision data which were recorded by the CDF II detector during Run II of 



S 



the Fermilab Tevatron. 



II. DETECTOR DESCRIPTION 

The CDF II detector is described in detail elsewhere {4]. The detector has a charged 
particle tracking system immersed in a 1.4 T solenoidal magnetic field coaxial with the 
p and p beams. Tracking over the radial range 1.5 cm to 28 cm is provided by three 
complementary silicon micro-strip detectors 5, [a, Q]. A 3.1-m-long open-cell drift chamber 
covers the radial range from 40 to 137 cm 8]. The fiducial region of the silicon system 
extends to pseudorapidity [9J] of \r]\ ~ 2, while the drift chamber (COT) provides coverage 
for 1 77 1 < 1. Segmented electromagnetic and hadronic calorimeters fiol. Q surround 
the tracking system and measure the energy of interacting particles in the pseudorapidity 
range \r]\ < 3.6. A set of drift chambers and scintillation counters [13] located outside the 
hadron calorimeter and another set behind a 60 cm iron shield detect muon candidates 
with \r]\ < 0.6. Additional chambers and counters detect muon candidates in the region 
0.6 < \rj\ < 1.0. Cherenkov counters [ljjl located in the 3.7 < \rj\ < 4.7 region measure the 
average number of inelastic pp collisions per bunch crossing and thereby determine the beam 
luminosity. 



III. DATA SAMPLE AND EVENT SELECTION 



In Run II of the Tevatron, protons and antiprotons collide at a center-of-mass energy of 
1.96 TeV. In such collisions, the SM predicts that ~ 85% of ti pairs are produced through 
quark antiquark annihilation and ~ 15% are produced by gluon fusion. Top quarks are 
expected to decay almost exclusively to a If boson and a b quark. The W subsequently 
decays either hadronically, to a pair of quarks or leptonically, to a lepton and a neutrino. The 
final state bblvqq' (where I = e, p only) resulting from one of each type of W decay is called 
the "lepton + jets" mode. This channel has a large branching fraction with a good signal- 
to-background ratio; we use it to measure m t using the decay length technique. Lepton + 
jets ti events typically contain a high transverse momentum (p?) electron or muon, missing 
transverse energy {St) from the undetected neutrino, and four high transverse energy (Et) 
jets, two of which originate from b quarks. All methods employed in previous measurements 
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of mt require the presence of all four jets, since they attempt to fully reconstruct the decays 
of both top quarks. Sometimes one or more of the jets may not be reconstructed, making 
those events ineligible for traditional methods. The decay length technique, however, can 
be applied to such events, providing the only measurement of m t from these data. 

Results reported here are obtained from an analysis of data collected between March 
2002 and September 2005. The data sample corresponds to an integrated luminosity of 
~ 695 pb _1 . CDF II employs a three-tiered trigger system to sequentially reduce event rates 
from ~ 1.7 MHz to ~ 80 Hz. The data used in this analysis were recorded with an inclusive 
lepton trigger that requires an electron (muon) with Et >18 GeV (pt >18 GeV/c). 

From this inclusive lepton dataset we select events offline with an electron (muon) with 
E T > 20 GeV (pr > 20 GeV/c), E T >20 GeV, and at least 3 jets with \q\ < 2 and energy- 



corrected [15[ Et >15 GeV 16]. Electron candidates are required to have a well-measured 
track pointing at an energy deposit in the calorimeter. The energy signature must be isolated 
from other calorimeter activity, and must have shower profiles consistent with expectation. 
We select muon candidates by requiring that the hadronic and electromagnetic energy de- 
posited by the candidate in the calorimeter be consistent with that of a minimum ionizing 
particle. In addition, we match partially reconstructed tracks in the muon chambers with 
well-measured tracks reconstructed in the COT. Jets are clustered using a fixed-cone algo- 



rithm with a cone size AR = J (Ar/) 2 + (A</>) 2 = 0.4. Finally, in order to better distinguish 
ti events from background processes, we require at least one jet in the event to be identified 
as a b jet ( "tagged" ) by the reconstruction of a secondary vertex within that jet as described 
in Section [TVl We refer to the dataset selected above as the "tagged lepton + jets" sample. 

IV. MEASUREMENT OF TRANSVERSE DECAY LENGTH 

The primary (primevtx) and secondary vertex algorithms (secvtx) used are described 



in Ref. 17|] . primevtx reconstructs the primary event vertex with a precision of ~ 15 /im in 
the plane transverse to the beam for ti events, secvtx exploits the relatively long lifetime 
of b hadrons in top decays to reconstruct a secondary vertex significantly displaced from 
the primary interaction. Secondary vertex tagging operates on a per-jet basis, where only 
tracks associated with the jet are considered. We require that these tracks have at least three 
silicon hits attached to them, are within 2.0 cm from the primary vertex in the longitudinal 
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direction, and that the final track fits have y^/d.o.f. < 8.0. We select tracks contained 
inside a jet which are displaced with respect to the primary vertex if they have a large well- 
measured impact parameter with respect to that same vertex. The secvtx algorithm uses 
a two-pass approach to find secondary vertices from these selected tracks. In the first pass, 
it attempts to reconstruct a secondary vertex which includes at least three tracks. If the 
first pass is unsuccessful, a second pass is attempted which makes tighter track requirements 
and tries to reconstruct a two-track vertex. 

Once a secondary vertex is found, we calculate the two-dimensional decay length, L 2D , 
as the projection onto the jet axis, in the plane transverse to the beam, of the vector 
pointing from the primary vertex to the secondary vertex. The sign of L 2 d is given by 
the 4> difference between the jet axis and the secondary vertex vector (positive if less than 
90 degrees, negative if greater than 90 degrees). The secondary vertices corresponding to 
the decay of b and c hadrons have large positive L 2 d, while the secondary vertices from 
mismeasured tracks form a Gaussian distribution centered around L 2 d = with a width 
corresponding to the detector resolution. A jet is "positively tagged" if its transverse decay 
length divided by the uncertainty on that measurement [L 2 d / 'o l 2D ) ls greater than 7.5. 
Similarly, a jet is "negatively tagged" if L 2D /cr L2D < —7.5. The positively tagged jet sample 
is predominantly composed of heavy-flavor (b or c) jets while the negatively-tagged sample 
is mostly composed of light quark jets. 

This analysis requires an accurate simulation of L 2 £>. To check the accuracy of the CDF II 
simulation we examine heavy-flavor enriched data samples. We use dijet data samples 
recorded with a trigger requiring an 8 GeV electron (muon). The lepton in these data, 
which we call a "soft lepton," often comes from the semi-leptonic decay of a b or c quark 
such that the heavy-flavor content of these samples is enhanced relative to generic dijet 
data. We compare these to herwig generic dijet MC samples which have been pre- 
selected to contain a soft lepton in analogy with the trigger requirements imposed upon the 
observed data. To increase the bb purity of these samples, we require that the soft lepton 
that triggered the event be contained within a jet with Et > 9 GeV that is also positively 
tagged by secvtx. We also require the presence of another jet with E T > 9 GeV, at least 
2.0 radians away in azimuth, that is also positively tagged. Finally, in order to distinguish b 
quarks from c quarks, we require the invariant mass of the four- vectors forming the secondary 
vertex 19] of the tagged jets to be greater than 1.5 GeV/c 2 . With this essentially triple- 
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tagged (two secvtx plus one soft lepton) selection, a purity of ~ 99% bb is obtained [20 ]. 
For all events passing the selection criteria, we make a histogram of the L 2 d for all positive 
tags. These histograms are shown with observed data and MC overlaid in Fig. [TJ We observe 
that the simulation models the L 2 d distribution very well. We quantify this agreement by 
comparing (L 2 d) for both observed data and MC as follows, where the errors are statistical 
only: 

(L% a ) = 0.378 ± 0.002 cm, (L%?) = 0.381 ± 0.004 cm (2) 

From these, we compute a data/MC scale-factor which could be applied to the mean trans- 
verse decay length measured in the observed data of 0.992 ±0.012. This number is consistent 
with a scale factor of unity; we conclude that our simulation models the transverse decay 
length of b hadrons with sufficient accuracy and do not apply any correction to the observed 
data. This ratio encompasses many different possible sources of discrepancy between our ob- 
served data and our MC simulation including effects from detector resolution, fragmentation, 
and the relative proportions and lifetimes of the various b hadrons. As such, it is a compre- 
hensive data-driven means of quantifying any systematic uncertainties in the measurement 
of (Z/2d). We note that in order to apply this treatment to b hadrons in ti events, we rely 
on the assumption of the universality of 6-fragmentation, i. e. that the fragmentation of a b 
quark is independent of the process in which that b quark was produced. This assumption is 



predicted by the QCD factorization theorem 



and is supported by a significant body of 



experimental evidence [221 ] . We assign the 1.2% statistical uncertainty on the central value 
in the above calculation as a systematic uncertainty on the accuracy of our MC simulation. 

We perform several additional checks on the data/MC scale-factor. Since the average 
energy of b hadrons from top decays is higher than that of those used to compute the 
scale factor, we examine the ratio as a function of jet Et- We find the scale factor to be 
independent of jet energy within uncertainties. We also compute a data/MC ratio of the 
(Z/2d) for negatively tagged jets, thereby more directly checking the resolution modeling of 
the simulation. We again measure a scale factor of unity within uncertainties. 
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FIG. 1: Comparison of L2D of positive tags in events with two secvtx tags and an identified 
electron or muon in one of the two tagged jets. Observed data and dijet MC are compared for this 
essentially triple tagged sample. Both the lepton tag and the non-lepton tag are included. Both 
distributions are normalized to unit area. 

V. SAMPLE COMPOSITION 

The tagged lepton + jets sample selected as described in Section [TV] has an expected 
signal-to-background ratio of about 2.5:1. The dominant background is the production of 
W plus multijet events. These events enter the signal sample when one of the jets is a b 
jet or c jet, or a light quark jet that has been mistakenly tagged as containing a secondary 
vertex. We call the latter type of events "mistags." The other substantial background 
comes from collisions which do not produce a W boson, termed "non-jy" events. These 
events are typically QCD multi-jet events where one jet has been misidentified as a high-p^ 
lepton and mismeasured energies produce apparent Additionally, other processes such 
as WW, WZ, ZZ, Z — > tt, and single top contribute small amounts to the tagged lepton 
+ jet sample. The techniques used to calculate the expected contributions to this sample 



are detailed elsewhere 
summarized in Table [B 



17| . Estimated contributions to the tagged lepton + jets sample are 
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TABLE I: Estimated number of events from background sources and single top that contribute to 
the tagged lepton + jets sample. The number of events observed data is also presented. The excess 
above the total background plus single top is assumed to be ti events. Errors include statistical 
and systematic uncertainties. 



Source 


Number of Events 


Wbb 


27.9 ±6.2 


Wcc 


12.2 ± 3.2 


Wc 


6.9 ±1.6 


non-W 


12.5 ±2.6 


Mistags 


40.9 ±3.8 


WW, WZ, ZZ, Z — > tt 


5.7 ±1.0 


Total Background 


106.1 ± 10.5 


Single Top 


5.3 ±0.5 


Total Background + Single Top 


111.4 ±11.0 


Data 


375 



VI. EXPECTED L 2D DISTRIBUTIONS 

[ — I 

We generate herwig ti MC samples using the CTEQ5L parton distribution functions [231 ] 
followed by a detailed simulation of the CDF II detector. The CLEO QQ Monte Carlo 
simulation models the decays of b and c hadrons. We produce these samples with top quark 
masses ranging from 130-230 GeV/c 2 in 5 GeV/c 2 intervals. We subject these simulated 
events to the identical event selection as that required of the observed data. After selection, 
we construct histograms of the transverse decay lengths of all positive tags in order to 
obtain L 2 d distributions for each mass point. A similar process is performed for each of 
the backgrounds described above. We model the L 2D distributions for the Wbb, Wcc, and 



Wc backgrounds using alpgen [24j matrix element calculations which have been interfaced 
with herwig to simulate the hadronization process. To model the L 2 d distribution from 
the W + mistag background we construct a hybrid data/MC template. Mistags can either 
come from tracks which appear displaced due to limited impact-parameter resolution or 
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from tracks that originate from actual long-lived particles that are not b hadrons. We use 
tags from negatively tagged W + jets data, reflected about L 2 d = 0, to obtain a data-driven 
positive mistag shape to model the resolution contribution to the distribution. We combine 
this with ALPGEN interfaced to herwig MC simulations of W + multijets, which we rely 
on to model the contribution from long-lived particles such as K$ and A. We obtain the 
relative normalization for this combination from independent studies comparing positive 
and negative tags [l^]. For the purpose of this analysis, WW, WZ, ZZ, Z —> tt events are 
considered mistags, and the mistag L 2D distribution obtained above is used to model their 
small contribution. We obtain the non-W background L 2 £> distribution directly from our 
observed data. We select events with identical criteria to those for the signal sample, except 
the requirement that the lepton be calorimetrically isolated, where instead we explicitly 
require the lepton be non-isolated. For most top analyses, single top is a background to 
the pair-produced signal. With the decay length technique, however, this is not the case. 
Although the correlation is not quite as strong as for pair- produced top quarks 25J, the 
(L 2 d) from b hadrons from single top decays is also correlated with mj. We use pythia 26] 
MC to model the single-top L 2 d distribution as a function of top quark mass from 130-230 
GeV/c 2 . 

As a cross-check on the modeling of L 2 £> distributions for the various background pro- 
cesses, we examine the observed data in the background-dominated one- and two-jet events 
in the lepton + jet sample. The L 2 £> distribution of positive tags in selected events from the 
one and two jet bins is shown together with expected signal and background contributions 
in Fig. [21 Good agreement between MC and experimental data is observed; a Kolgomorov- 
Smirnov (KS) test yields a p-value of 30.6%. 



VII. ESTIMATION OF TOP MASS DEPENDENCE ON L 2D 

We treat the signal and background L 2D distributions described in the previous section 
as probability density functions from which we form ensembles of simulated experiments. In 
forming each ensemble, the number of events from a given background source (and single 
top) is obtained by allowing the number of events for each process to fluctuate separately 
about the expected contributions listed in Table [B The number of events from ti, which is 
similarly allowed to fluctuate, is taken to be the excess of the observed data in the tagged 
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FIG. 2: The L 2 d distribution of positive tags in lepton + jet events with 1 or 2 jets, for which only 
a small top contribution is expected. The points are the observed data. Expected contributions 
from signal and background MC are displayed cumulatively in the histogram. To facilitate the 
shape comparison, the MC is normalized to the observed data. 



lepton + jets sample over the summed contributions of the background processes and single 
top production. 

In computing (L^d), we will sum over tags rather than events. We convert the number 
of events for each process to a number of tags by multiplying by the probability, obtained 
from MC simulation, for that process to contain more than one SECVTX tags. We repeat 
this procedure 1,000 times for each mass point over the full mass range 130-230 GeV/c 2 . We 
construct histograms of the {L 2 d) that results from each pseudo-experiment performed at 
a given m t . We extract the mean and ±lcr variance from these histograms, for each value 
of mt, and fit these points to third-degree polynomials. The fit to the mean establishes the 
most probable value for a true top mass given a measured (L 2 d) and is the function used 
to make the top mass measurement from the (L 2 d) observed in the data. Similarly, the fits 
to the variance form ±lcr Neyman 27j confidence intervals which will be used to give the 
statistical uncertainty of the measurement in Section HXl as shown in Fig. [31 
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We derived these functions and validated our method prior to examining our experimental 
data. We employed simulated data ensembles containing MC it events with unknown top 
quark mass to demonstrate that m t could be extracted accurately and with appropriate 
precision. 



VIII. ESTIMATES OF SYSTEMATIC UNCERTAINTIES 



The systematic uncertainties for this measurement come from three kinds of sources. The 
first arises from the accuracy of the modeling of factors which affect the top (or subsequent 
bottom) quark's momentum. We estimate the uncertainty due to initial- and final-state 
gluon radiation by varying the relevant parameters by ±1 standard deviation in the simula- 
tion [28] and observing the effect on the measured m t . We quote half the difference between 
these variations as the systematic uncertainty. To assess the systematic uncertainty from our 
choice of the cteq5l parton distribution function, we observe the shift in m t that results 
when we substitute the mrst72 and mrst75 sets [29] which are evaluated at different values 
of the str ong coupling constant, a s . Additionally, we vary the 20 eigenvectors in the cteq6m 



package 



30] by ±1 standard deviation. We add the shifts observed in mt in quadrature to 



determine the total systematic uncertainty. We estimate the systematic uncertainty due to 
our choice of herwig to simulate our signal events by measuring the shift induced in m t 
upon substitution of pythia simulations. Finally, we evaluate the uncertainty due to the 
jet energy scale by applying jet energy corrections which have been shifted by ±1 standard 
deviations [l6] and taking half the difference as the systematic error. As anticipated, we 
find the analysis to have negligible sensitivity to such shifts. 

The second type of systematic uncertainty comes from potential inaccuracies in the size 
or shape of background L 2 d distributions. We quantify the uncertainty due to background 
normalization by increasing/decreasing the contributions from each background process ac- 
cording to the uncertainties listed in Table [B We estimate the effects of uncertainties in the 
shape of the background L 2 d distributions by substituting altered distributions and noting 
the corresponding shift in m t . For the W + heavy-flavor background shapes we vary the 
momentum transfer parameter (q 2 ) in the calculation. For the non-W shape we substitute a 
distribution obtained from observed bb data. Finally, to estimate the systematic uncertainty 
due to the distribution used to model the mistag background, we measure the effect on m t 
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TABLE II: Summary of sources of systematic error and their estimated uncertainties. 



Source of Systematic Error 


Uncertainty (GeV/c 2 ) 


Monte Carlo Generator 


0.7 


Initial State Gluon Radiation 


1.0 


Final State Gluon Radiation 


0.9 


Parton Distribution Functions 


0.5 


Event Selection (Jet Energy Scale) 


0.3 


Background Shape 


6.8 


Background Normalization 


2.3 


Multiple Interactions 


0.2 


Data/MC (L 2D ) Ratio 


4.2 


Total 


8.6 



when we alternatively use, only the shape derived from the observed data, and only the shape 
derived from MC, in place of the data/MC hybrid distribution used in the analysis. We take 
half the difference between these two determinations as the systematic uncertainty. We add 
the separate background shape uncertainties in quadrature to arrive at a total systematic 
uncertainty. 

The final type of systematic uncertainty comes from imperfections of detector simulation 
of L 2 d or other, experimentally indistinguishable, disagreements between the (L 2 d) in MC 
and observed data that may arise from inaccuracies in the simulation of b hadron decays. 
This uncertainty is taken as the error on the (L 2 d) data/MC ratio as discussed in Section HVl 
Table [II] lists each of the sources of systematic error and their corresponding uncertainties. 
A total systematic uncertainty of 8.6 GeV/c 2 is assigned to the measurement. The dominant 
sources of systematic error are the finite statistics used to derive the (L 2 d) data/MC ratio 
and the modeling of the L 2 d distribution from mistags. 
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IX. RESULTS 



From 456 positive secvtx tags in 375 events in the lepton + jets sample corresponding 
to 695 pb _1 we measure 

(L 2D ) = 0.581 cm (3) 

We draw a vertical line at this (L 2 d) and read off the intersections with the most probable 
value and ±lcr confidence interval curves obtained above as illustrated in Fig. [3] to extract 
a measurement of 

m t = 180.7 tm ( stat ) ± 8.6 (syst.) GeV/c 2 . (4) 

The L 2 d distribution of positive tags in selected events, from which the mean transverse 
decay length used to measure top mass is extracted, is shown together with expected contri- 
butions from signal and background MC overlaid in Fig. HI Reasonable agreement between 
data and MC is observed; a KS test yields a p-value of 16.7%. 




0.48 0.5 0.52 0.54 0.56 0.58 0.6 0.62 0.64 

<L 2D > [cm] 



FIG. 3: Most-probable (solid) and ±l<r (broken) mt curves as a function of mean transverse decay 
length. Uncertainties are statistical only. Measured mean transverse decay length is overlaid as 
dashed line. 
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FIG. 4: The L2D distribution of positive tags in selected events from which the mean transverse 
decay length used to measure top mass is extracted. The points are the observed data. Expected 
contributions from signal and background MC are displayed cumulatively in the histogram. To 
facilitate the shape comparison, the MC is normalized to the observed data. 



X. CONCLUSION 



We have performed the first measurement of the top-quark mass using the decay 
length technique. Using 695 pb -1 data from Run II of the Tevatron we measure m t = 
180.7 -13.4 (stat.) ± 8.6 (syst.) GeV/c 2 , consistent with the SM expectation as well as 
the combined result from of all existing measurements, m t = 171 A ± 2.1 GeV/c 2 (31 ]. 
Since it has negligible dependence on the jet energy scale and analyzes lepton plus three jet 
events which are not used by any other top mass measurement, the decay length technique 
is largely uncorrelated [32j with other methods. Consequently, while the result presented in 
this paper is not a competitive measurement of m t by itself, it will help to reduce the overall 
uncertainty on the top mass when combined with other results. In the combination refer- 
enced above, this measurement contributes at the few percent level to the overall combined 
result. 

This measurement is statistically limited and its dominant systematic uncertainties are 
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likely reducible. The precision of this measurement, therefore, will continue to improve over 
the course of Run II of the Tevatron. The asymptotic performance of this technique at the 
Tevatron (and the LHC) are studied in detail in [2}. The current measurement establishes 
the technique and supports that paper's claim that the decay length method is a useful 
complement to existing measurements at the Tevatron and beyond. 
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